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Abstract — For  interferometric  SAR  (InSAR)  processing,  the 
precision  of  the  interferometric  phase  is  very  important  since  it 
decides  directly  that  of  the  final  measurement.  It  is  known  that 
the  interferometric  phase  is  composed  of  many  different  parts 
such  as  the  phase  due  to  the  topography,  the  possible 
deformation,  the  atmosphere,  the  track  error  and  so  on.  Before 
doing  phase-unwrapping,  the  interferogram  must  be  flattened  to 
decrease  the  density  of  the  fringes  in  the  interferogram  or  it  may 
make  phase  unwrapping  fail.  For  the  practical  application  of 
InSAR,  if  the  interferogram  is  not  flattened  well,  the  residual 
phase  due  to  the  earth  will  be  considered  to  be  the  phase  due  to 
the  topography  or  the  deformation.  In  this  paper  the  different 
earth  reference  models  that  have  been  used  for  doing  flattening 
are  analyzed.  Some  experiments  have  been  done  on  the  real 
images  and  the  conclusions  are  given  at  the  end. 

I.  Introduction 

The  application  of  the  InSAR  technique  is  greatly  related  with 
the  accuracy  of  the  interferometric  phase  which  is  decided  by 
the  coherence  of  the  images  data.  When  a  pair  of  SAR  images 
is  well  registered,  their  phase  difference  images  called 
interferometric  phase  or  interferogram  will  show  the  fringes 
which  are  related  with  the  information  of  the  topography  on  the 
earth.  It  is  known  that  the  interferometric  phase  is  composed  of 
many  different  parts  such  as  the  phase  due  to  the  topography 

(f>topo  •>  the  flat  earth  (f)earth  ,  the  possible  deformation  (j)dej. ,  the 
atmosphere  (j)atm  ,  the  track  error  (f>track  and  so  on  [2]. 


t  -  topo  +  dearth  +  tatm  +  tdef  +  track  +  t noise  ( 1 ) 

If  it  is  not  done  properly,  the  phase  due  to  the  manually 
processing  will  be  introduced  into  the  final  inteferometric 
phase.  It  has  been  shown  that  the  precision  of  the 
interferometric  phase  decides  directly  that  of  the  final 
measurements.  Though  the  quality  of  the  interferogram  suffers 
mainly  from  the  coherence  of  the  data,  the  method  selected  for 
interferogram  processing  is  also  very  important.  A  good 
interferogram  generation  method  will  greatly  reduce  the  effect 
due  to  the  manually  processing. 

According  to  the  InSAR  imaging  geometry,  the  raw 
interferogram  can’t  avoid  the  flat  earth  effect  that  means  the 
phase  difference  created  due  to  the  non-zero  baseline.  Because 
of  the  existence  of  the  flat  earth  effect,  the  objects  with  the 
same  height  on  the  ground  don’t  have  the  same  interferometric 


phase.  Furthermore  the  flat  earth  effect  makes  the  density  of 
the  fringes  in  the  interferogram  very  high,  and  it  adds  more 
difficulty  for  the  vital  phase-unwrapping  processing.  So  before 
doing  phase-unwrapping,  the  interferograms  must  be  flattened. 
At  the  mean  time,  after  doing  flattening  the  interferogram  will 
look  like  the  contour  image  of  the  topography. 

II.  Problem  Analysis 

The  flat-earth  removal  methods  can  be  categorized  into  two 
classes:  using  the  parameters  of  the  orbit  and  dealing  with  the 
interferograms  directly.  The  problem  with  the  first  method  is 
the  uncertainty  of  the  orbit  parameters  will  have  an  effect  on 
the  positioning  accuracy.  The  second  method  is  mainly 
realized  by  estimating  the  position  of  the  maximum  frequency 
in  the  interferogram  spectrum  where  the  regular  flat  earth 
phase  will  dominate.  Then  the  reference  phase  can  be 
removed.  But  this  kind  of  method  suffers  from  the  noise  due 
to  the  low  coherence  of  the  data.  For  those  low  coherent  SAR 
images  data  especially  those  that  have  long  time  baseline,  this 
kind  of  problem  is  very  seriously.  The  flattening  results 
depend  greatly  on  the  maximum  estimation  accuracy. 

The  orbit  parameters  based  interferogram  flattening  methods 
are  mainly  using  the  flat  earth  model  or  the  elliptical  earth 
model. 


Figure  1 .  Flat  earth  based  flattening  model 


Fig-1  illustrates  the  flattening  method  using  the  flat  earth 
model,  where  B  is  the  baseline,  6  is  the  incidence  angel, 
OC  is  the  tilt  angel,  \\  and  r2  are  the  distance  between  the 

satellite  and  the  object  on  the  earth,  H  is  the  height  of  the 
satellite  and  h  is  the  height  of  the  object  on  the  earth.  The  flat 
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earth  phase  can  be  calculated  using  (2). 


(6) 


dearth  =  "^  ^  5  SUl(0  -  d)  (2) 

While  for  the  flat  earth  model,  it  is  easy  to  know  sin  6  and 
COS  0  by  equ  (3) 


c„s«  =  ^ 


Then  (pearth  can  be  rewritten  as: 


^  earth 


4  K 

T 


B 


H2  .  H 
coso.  1 - - — sin  a — 


V 


(3) 


(4) 


For  the  elliptical  earth  based  flattening  method,  before  doing 
flattening,  the  position  of  each  pixel  in  the  image  has  to  be 
calculated  using  the  method  proposed  by  [1],  Because  the 
satellite  coordinate  can  be  calculated  using  the  orbit  data,  it  is 

easy  to  calculate  the  flat  earth  or  reference  phase  (j>earth . 

4  Jl 

dearth  (5) 


AO  = - ,  h  «  R,  &  R 

R{  sin  0 

Let  A  =  6  —  OC  ,  then 

2?sin(A  +  A#)  =  B  sin  A  cos  A&  +  B  cos  A  sin  AO  (7) 


Because  COS  AO  ~  1  and  sin  AO  ~  AO  ,  equ  (7)  is  equal  to 
equ  (8) 

h 

5sin(A  +  A#)  =  i?sin  A  +  5cos  A -  (8) 

Rx  sin  0 


From  (8),  the  phase  difference  (j)^  due  to  the  application  of 
different  earth  reference  model  is: 


hw  =^-Bcos(0-a) 


R  sin  0 


(9) 


The  phase  difference  (f)^  can  be  related  with  the  registration 
error  pdiy  by  (10),  where  Rs  is  the  slant  range  resolution. 

h 


Pdiff  =  Bcos(O-a) 


RxRs  sin  6 


(10) 


The  flat  earth  effect  can  be  removed  point  wisely.  But  it  will 
cost  much  more  time  when  the  image  size  is  very  large.  In 
order  to  save  time  a  polynomial  can  be  used  to  fit  the  reference 

phase  (j)earth  •  It  should  be  a  function  of  the  pixel  coordinate 

(row,  column),  which  means  it  changes  both  in  range  direction 
and  in  azimuth  direction. 

From  fig-3  it  is  easy  to  get  (6). 


CONCLUSIONS 

For  the  different  earth  reference  models  that  have  used  to  do 
interferogram  flattening,  the  computation  load  is  a  little 
different.  While  the  selection  of  the  reference  earth  model  will 
depend  on  the  flying  height,  the  size  of  the  imaging  area  and 
the  incidence  angel  of  the  synthetic  aperture  radar.  The  flat 
earth  model  based  flattening  model  will  run  the  risk  of 
decreasing  the  SAR  images  fine  registration  precision  which 
will  have  an  impact  on  the  coherence  of  the  data. 
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Figure  4.  Flattening  results  using  Beijing  area  ERS— 2  May  20th  ,1998  and  Jun  24th  ,  1998  (bottom  :  amplitude  image;  middele:  polynomial  fitted  reference 

phase;  top:  flattening  result  based  on  elliptical  model) 


